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Abstract 

After discussing results of dilepton rates and electrical conductivity obtained from continuum extrapolated results 
of light quark correlation functions in quenched QCD I will give a review on recent developments in lattice QCD 
studies for quarkonia in the quark gluon plasma. Recent progress in the extraction of spectral properties from lattice 
QCD calculations of hadronic correlation functions will be discussed. Besides medium modifications of bound states 
and their dissociation in the plasma I will focus on transport coefficients, like heavy quark diffusion constants extracted 
from different correlation functions on the lattice. Present limitations and future perspectives for studies of quarkonia 
and related transport coefficients on the lattice will be discussed. See also 0] El El [4 1 for recent overviews and related 
discussions of those topics. 
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1. Dilepton rates and electrical conductivity 

A detailed knowledge of the spectral function of mesonic states at high temperatures is of fundamental importance 
for the understanding of the properties of the quark gluon plasma as well as for the analysis of this hot and dense 
medium in heavy ion collision experiments. In lattice QCD calculations the spectral function itself is not directly 
calculable but can be derived from hadronic correlation functions. In this section we will focus on the vector channel in 
the light quark sector before discussing heavy quark results in the remaining sections. The current-current correlation 
functions can be represented in terms of an integral over spectral functions, p^fa),^, T). We denote by p„ the sum 
over the three space-space components of the vector spectral function and also introduce the vector spectral function 
p v = poo + pa. With this we obtain the corresponding correlation functions, 
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du> . cosh(oXr - 1 /2T)) 

Gh(t, P,T)= — p H (co, p, T) . \ ' , H = 00, ii, V. (1) 

1 2n sinh(w/2r) 



The vector spectral function is directly related to the thermal production rate of dilepton pairs with squared invariant 

mass to 2 - p 2 , 

_ r aL p v (u),p,T) 

dojtfp Pm 67T 3 (CO 2 - p 2 )^ 1 - 1) ' ' 

where a em is the electromagnetic fine structure constant. This expression is correct to leading order in the electromag- 
netic interaction but contains the full QCD information once the vector spectral function is known, (see also Q). In 
the following only the p — case will be discussed. 
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Figure 1: Lattice QCD results together with the continuum extrapolation of the vector correlation function at 1.46 T c from | 5 1. 



In a controlled continuum extrapolation of the vector correlation function calculated in quenched QCD at 
T - 1.45 T c was performed for the first time. Those lattice results of |5| together with the result in the continuum 
limit are shown in Fig.[T] 

The vector correlation function in the continuum was used to fit to a phenomenological inspired Ansatz for spectral 
function, 

poo(w) = -2nx q u>6(tj) , (3) 

Pii (u) = 2 Xq c BW —^^—r + 2- (1 + k) co 2 tanh(w/4r) , (4) 
oj l + (r/2) z 2n 

which contains a Breit-Wigner contribution at small frequencies that gives the correct w-dependence in the small 
frequency limit to determine the electrical conductivity and a correction parameter k(T) that parameterizes deviations 
from a free spectral function at large energies. The result of this fit shown in Fig. [T] demonstrates that already this 
minimal ansatz provides a good description of current numerical results for the vector correlation function. 




Figure 2: Thermal dilepton rate in 2-flavor QCD (left). The HTL curve is for a thermal quark mass mjjT = 1 and the Born rate is obtained by 
using the free spectral function. The right hand part of the figure shows the spectral functions that entered the calculation of the dilepton rate, (see 
text and 1 5 1 for details). 

Using Eq. |2} the continuum estimate for the dilepton rate at T - 1.45 T c in Fig.|2jleft) are obtained. They merge 
with the Born-rate at a> > 4 and are well above this rate for lower energies, still staying below the HTL result at 
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Figure 3: Pseudo-scalar screening mass (a) and amplitude (b) divided by the corresponding zero temperature values as function of the temperature. 
The solid black line is the free quark result. The plots are from (9). 



small a). In contrast to the HTL result, the limit a> — » is well behaved (see Fig. fright)). In this limit the spatial 
components of the vector spectral function determine the electrical conductivity 

cr C em pa(u>) 

— = hm . (5) 

T 6 w-»0 uT 

The vector spectral function shown in Fig.|2]at vanishing energy leads to an estimate of the electrical conductivity 

1/3 < —^—^ < 1 at r-1.45T c . (6) 

Although only a lower and upper limit was estimated, the result demonstrates that transport coefficients can be reliably 
extracted from continuum extrapolated lattice QCD correlation functions which in this respect improves on earlier 
results J6l|7l. See J5 1 f° r more details and [8 1 for a recent study using a new technique to extract the spectral function 
in this limit. 

The vector spectral function at light-like 4-momentum yields the photon emission rate of a thermal medium, 

d^y _ r a em p v (<jj -\p\,T) 
"dip' em 4*2 eflT-l ■ U) 

The emission rate of soft photons, thus can be related to the electrical conductivity, 

dR y 3 

hm cd— = —cr(T)Ta em . (8) 
w-»o Gfp 2n l 

Together with Eq. (|7| this yields for the zero energy limit of the thermal photon rate, 

lim = (0.0004 - 0.0013) T 2 C * (1 - 3) • 10~ 5 GeV 2 at T * 1.45 T c . (9) 

w^o d j p 

Results of the vector correlation function at non-zero momentum will be used in the future to estimate p + dilepton 
rates and the finite energy contribution of photons in the QGP. 

2. Charmonium screening masses 

High precision data at a large number of Euclidean time separations are needed to extract spectral properties from 
temporal correlation function. An alternative approach that may provide additional information on the in-medium 
modifications of spectral functions and medium effects of mesonic states are spatial correlation functions, 

G(z,T) = Jdxdyj^ dT{J(x,y,z,T)J(0, 0,0,0)), (10) 
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Figure 4: Results for the charmonium spectral functions in the vector channel. Shaded areas are statistical errors from Jackknife analyses. 
Horizontal error bars at the first peaks of spectral functions at 0.73 T c an 1.4 T c indicate the statistical uncertainties of the peak location obtained 
from the Jackknife analyses, (see 1 1 1 for details). 



which are related to the momentum dependent spectral function through 

Jo 2k lo J_ M 2n 



(11) 



At large distance, zT > 1, where medium effects are expected to be largest, G(z, T) decays exponentially allowing to 
define a screening mass M scr . If bound states appear in the spectral function, the screening mass is determined by the 
mass M of the lowest lying meson state. On the other hand for unbound quark antiquark pairs at high temperature, 
M scr is given by the lowest non-vanishing Matsubara frequency 2 ^J{nT) 2 + m\ where m c is the quark mass. Those 
two limits and the transition between them may serve a an indication for significant modifications and dissociation of 
heavy quark meson states. 

In Fig.[3]the results for the charmonium screening masses calculated for a dynamical 2+1 flavor lattice QCD from 
||9l are shown, normalized to the zero temperature value. The change in the temperature dependence becomes apparent 
around 1 .5 T c where deviations from the zero temperature limit become large and the behavior becomes compatible 
with that of unbound cc pairs. This is an indication that charmonium states melt around this temperature which is a 
consistent picture with the findings from the study of temporal charmonium correlators in quenched QCD [10] that 
will be discussed in the following section. 



3. Charmonium spectral functions and heavy quark diffusion 

The properties of charmonium states at finite temperature were studied in [10] in quenched lattice QCD on large 
and fine isotropic lattices. Although no continuum extrapolation was performed, the results for lattice sizes of 128 3 xJV T 
for N T between 16 and 96 and lattice spacing a down to 0.01 fm are already close to the continuum and cut-off effects 
in the spectral function are well separated from the physically interesting energy regions, as am c is much smaller 
than one. In IfTOll a detailed analysis of statistical and systematic uncertainties and default model dependencies in 
the maximum entropy method (MEM) used to determine the spectral functions from the lattice calculations of the 
charmonium correlation functions was performed. 

The results of this study, shown in Fig. [4] for the vector channel, suggest that charmonium states are dissociated 
already at 1.46 T c . While at 0.73 T c stable and reliable ground state peaks are found, at the higher temperatures in this 
study no signals for bound state contributions but rather a threshold enhancement are seen. 

The slope of the vector spectral function in the w — > limit again determines a transport coefficient, the heavy 
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Figure 5: Spectral functions p(to) in the vector channel (T), normalized with the heavy quark mass. This plot is from II 11 . 



quark diffusion constant D, 

1 v-i pUoi,P=0,T) 
D = —r, lim > " — , (12) 

/L i= 1 

where x°° is the quark number susceptibility. While compatible with zero below T c , almost temperature independent 
values of 2nTD around two at T — 1.46, 2.20 and 2.93 were obtained in IfTUl (see Fig.|7J. 

This is a further indication of strong medium effects and Jf¥ being melted at 1.46 T c and higher temperatures. 
Together with the results for the charmonium screening masses of the previous section this gives a consistent picture 
and shows that the interesting temperature region for medium effects and the study for dissociation of charmonium is 
between T c and 1.5 T c . Future studies at temperatures closer to T c using continuum extrapolated correlation functions 
may allow a more detailed determination of the dissociation patterns of the different charmonium states in the quark 
gluon plasma. 



4. Bottomonium from lattice NRQCD 

At the energies probed at the LHC bottomonium spectroscopy becomes an important tool for the analysis of the 
produced hot medium. It may offer a cleaner probe compared to charmonium as feed-down effects from heavier 
quarkonium states are missing and regeneration from the thermal medium will be negligible. Furthermore the disso- 
ciation temperatures are far higher than for charmonium and therefore well separated from the hadronization stage. 
First results from the CMS experiment at LHC show a suppression of the higher excited states T{2S +3S) lfT2l and 
already indicate the importance of bottomonium in this respect. Therefore a detailed knowledge on the behavior and 
dissociation temperatures of bottomonium states is of fundamental importance. 

For bottom quarks the heavy quark mass scale can be integrated out in terms of an effective field theory, non- 
relativistic QCD (NRQCD). In a lattice discretization of this effective theory where the physics above the scale mt, 
is encoded in the perturbatively determined coefficients of the NRQCD Lagrangian, the spatial lattice spacing acts 
as a short distance cut-off and must satisfy mba s > 1 iTTTl IT31 IT4ll . The drawback of this is that no continuum limit 
is possible in lattice NRQCD and furthermore that cut-off effects are present in the spectral function already in the 
energy region relevant for bottomonium excited states as the cut-off effects and also the accessible energy region is 
determined by the spatial lattice spacing. 
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Figure 6: Left: Comparison of lattice with the NLO weak-coupling expansion [15]. Right: Like on the left but normalized to the model of eq. |15| 
Closed symbols correspond to k = 4T 3 , open ones to T = 2.25T C (see 1 16 1 for details). 



In contrast to a computationally expensive inversion of the fermion matrix, NRQCD propagators are calculated 
as an initial value problem. Information on transport properties are exponentially suppressed which is a further 
consequence of the quark mass scale acting as a low energy cut-off at 2mj. Nevertheless this facilitates the extraction 
of spectral properties in the bound state region as no constant contribution appear in the NRQCD correlation functions. 

Despite the limitations of the method, the results of IfTTTl demonstrate the benefits of lattice NRQCD. Their results 
for the spectral function in the vector channel (T) are shown in Fig. [5] It suggests that the ground state survives up to 
the highest temperature analyzed in this work and indicates the melting of higher excited states at temperatures above 
1 .4 T c , See [11] for a detailed discussion of the results. 

Especially for the behavior of the excited states a detailed analysis of the systematic uncertainties of this approach 
is still needed, especially to understand the cut-off effects and default model dependencies of the maximum entropy 
method used to extract the spectral functions. In addition lattice NRQCD results will be useful for potential model 
studies of quarkonium systems and will give additional input for future full relativistic lattice QCD calculations of 
bottomonium states. 



5. Heavy quark momentum diffusion coefficient 

Effective field theory methods are also important to derive operators to extract transport coefficients from lattice 
QCD. In the large quark mass limit Heavy Quark Effective Theory (HQET) can be used to derive a purely gluonic 
operator which is related to the "momentum diffusion coefficient" in the low-frequency limit, 

2Tp E (cu) 

k = hm . (13) 

fc)-»o a> 

Heavy quarks carry a color charge and are therefore subject to a colored Lorentz force. Like with other transport 
coefficients the corresponding "low-energy constants" are easiest to define at vanishing three-momentum when the 
Lorentz-force is proportional to the electric field strength. This leads to a "color-electric correlator" lfT71fT8l . 



U(};T)gE i (T,0)U(T;0)gE i (0,Q) 



1 3 (ReTr 

GfE( T ) = --y 

J i=i (ReTr[£/(i;0)] 



(14) 



where gE, denotes the color-electric field, T the temperature, and U(j2\ t\) a Wilson line in Euclidean time direction. 
If the corresponding spectral function, can be extracted lfl9l . then the "momentum diffusion coefficient", often 
denoted by k, can be obtained from Eq. ( p"3j ). According to non-relativistic linear response relations (valid for M » nT, 
where M stands for a heavy quark pole mass) the corresponding "diffusion coefficient" is given by D = 2T 2 /k. 
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Figure 7: Comparison of various estimates for the heavy quark diffusion constant (see the text for more details). 



Due to the gluonic nature of the operator Eq. ( p"4| ) large fluctuations lead to weak signal/noise ratio. Using clever 
improvement algorithms in addition with tree-level improvement to reduce cut-off effects and a perturbative renor- 
malization of the operator leads to results that show almost no lattice artifacts. 

In Fig. |6jleft) the results of |[T6l are shown in comparison with a NLO weak-coupling expansion, both normalized by 
the leading order [ 1 8 1 correlation function. To extract the diffusion coefficient we use a model for the spectral function 

I cjk\ 

Pmodei(w) = maxjp NL0 (w), — > , (15) 

with the free parameter k representing directly the momentum diffusion coefficient according to Eq. (fl3j, and the 
corresponding Euclidean correlator computed from 



7 model \ 



~ Jo * 



du , A Cos h((i-rr )f) 
sinh(#) 



(T) = Pmodel(w) ; \ ■ (16) 



Results for the ratio of the lattice data and the model correlation function for different values of k are shown in Fig. [6] 
and allow for an estimate for the momentum diffusion constant at the three temperatures analyzed. Converting those 
estimates to the usual diffusion coefficient D = 2T 2 /k one obtains DT ~ (0.5 ■ • ■ 0.8) which is slightly larger than the 
values obtained for the charm case [10] (see section [3} and on this qualitative level consistent with a similar analysis 
in (2D). 

Although the values for k and DT are to be understood as a first rough estimate, it indicates that a consistent picture 
appears to emerge. Improving the lattice data for the correlation function and performing the continuum limit, which 
is possible in this approach, will allow to test improved models for the spectral function and use methods like the 
maximum entropy method or other sophisticated methods to perform the analytic continuation like the one proposed 
and tested in I12T1I81. 



6. Conclusions and outlook 

Lattice QCD calculations have proven to provide important information on the in-medium modifications of hadronic 
states for light up to bottom quarks. High quality results and sophisticated methods to determine spectral properties 
from lattice correlation functions allow a reliable estimate of transport coefficients. This is accompanied by effective 
field theory methods used in the heavy quark mass limit. Especially for the analysis of bottomonium and even more 



O. Kaczmarek / Nuclear Physics A 00 (2012) 743] 



8 



for the determination of transport coefficients from lattice correlation functions additional information can be gained 
from this combination. 

Although most results are still obtained in the quenched limit or for rather large light quark masses, a consistent 
picture of medium modifications and the dissociation of quarkonia as well as of transport coefficients needed to 
describe the thermal medium and the behavior of heavy quarks in this medium emerges. 

For future charmonium studies the interesting temperature region between T c and 1.5 T c shall be explored in 
the future as the recent results indicate that charmonium is melted already above those temperatures. At least for 
the bottomonium ground states higher temperatures seem to be relevant and full relativistic lattice calculation are 
necessary to confirm the findings of the exploratory NRQCD study in this sector. 

In Fig. [7] we compare the values for the diffusion coefficients obtained from the analysis of the charmonium 
spectral function in quenched QCD IfTOll and the heavy quark momentum diffusion coefficient in the large quark mass 
limit [16 1 (light blue band). On this qualitative level all results are compatible with estimates derived from the heavy 
flavor Raa and v w PHENIX data El l23l . See also [24] for a determination of the diffusion constant using heavy 
quark free energies in a T-matrix approach and [25 1 for a recent discussion of heavy quark diffusion constants. Just 
for comparison we also show the result obtained from the electrical conductivity in the light quark sector Il5l l26l . 

The results of section[4]and[5]show the opportunity and applicability of effective field theory methods and operators 
derived from it for the study of heavy quark spectral and transport properties in lattice QCD calculations. Despite their 
limitations they will allow to constrain models, e.g. for potential model calculations J3]|27), and contribute additional 
information and input needed in full relativistic lattice QCD calculations in the future especially for the difficult task 
of extracting spectral properties from hadronic correlation functions. 

In a similar way as for the light quark case discussed in section [T[ an extrapolation of the correlation functions to 
the continuum will be required in the heavy quark sector to remove any lattice discretization effects. Although at the 
moment only possible in quenched QCD this will allow a more quantitative determination of spectral and transport 
properties for quarkonium systems. This will lay the foundation for the analysis of quarkonium correlation functions 
in full QCD with realistic light quark mass becoming available in the coming years at least for moderate temporal 
lattice extensions and will allow to determine the effects of light quark masses on spectral properties of heavy mesons. 
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